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1.  INTRODUCTION 


An  optical  background  suppression  technique  is  described  which  greatly 
enhances  the  detectability  of  localized  sources  in  the  presence  of  a  nonuniform 
field  of  background  radiation.  The  technique  allows  both  the  spatial  and 
spectral  characteristics  of  the  source  to  be  exploited  in  order  to  aid  the 
detection  process,  while  it  suppresses  the  contribution  of  the  background  to 
the  observed  signal.  The  optical  system  consists  of  an  imaging  double-beam 
i nterferometer-spectrometer  which  has  been  configured  to  operate  also  as  a 
spatial  filter  by  suitably  tailoring  the  modulation  transfer  functions  (MTF) 
of  the  two  main  optical  paths  through  the  interferometer.  In  its  role  as  a 
spatial  filter,  the  double-beam  interferometer  preferentially  suppresses  the 
lower  spatial  frequencies  found  in  the  scene.  The  type  of  background  suppression 
discussed  here  is  especially  important  in  the  remote  detection  of  small  stationary 
sources  against  a  strong  spatially  structured  background  in  the  infrared  region 
of  the  spectrum.  In  particular,  the  technique  minimizes  the  image  processing 
required  by  significantly  reducing  the  background  before  the  image  is  recorded 
electronically. 

The  theory  of  the  imaging  double-beam  interferometer-spectrometer  with 
beams  having  suitably  matched  (tailored)  modulation  transfer  functions  (MTF)  is 
presented  in  §  2.  A  number  of  systems  with  tailored  MTFs  are  described,  and 
their  response  to  a  variety  of  spatially  structured  backgrounds  are  calculated. 

A  laboratory  double-beam  interferometer-spectrometer  with  tailored  MTF 
(defocussed  system)  is  described  in  §  3.  Details  of  the  optical  system,  mirror 
drive,  spectrum  analyzer,  and  calibration  facility  are  presented. 
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A  series  of  laboratory  tests  and  demonstrations  with  the  double-beam 
interferometer-spectrometer  have  been  carried  out.  These  tests  are  discussed 
in  4  and  include  the  response  of  the  interferometer-spectrometer  to  a  range 
of  spatial  frequencies  and  the  use  of  high  resolution  spectral  discrimination 
in  suppressing  backgrounds. 

Several  applications  of  the  optical  background  suppression  technique  are 
proposed  in  §  5.  These  include  the  remote  sensing  of  pollutants  in  smokestack 
effluents,  and  the  remote  sensing  of  missile  trails  from  surveillance  satellites 
The  various  sections  of  this  report  have  been  prepared  so  that  each  may 
be  read  independently  without  detailed  study  of  the  other  sections. 


2. 


PRINCIPLES  OF  THE  BACKGROUND  OPTICAL  SUPPRESSION  SCHEME 


2 . 1  Double-Beam  In  ter f e rometers 

The  technique  of  double-beam  interferometry  was  first  suqqested  by  Fellgett1 
in  1957.  Subsequently,  double-beam  interferometers  in  various  forms  have  been 
used  for  laboratory  and  astrophysical  measurements  in  which  the  spectra  to  be 
measured  were  superimposed  on  large  background  signals''-7.  In  these  measure¬ 
ments,  essentially  one  input  beam  contained  radiation  from  the  source  and  back¬ 
ground  while  the  other  input  beam  contained  only  background  radiation.  This 
has  been  done  by  using  adjacent  fields  of  view  with  a  nearly  uniform  background, 
by  simulating  a  well  known  background  (e.g.,  blackbody),  or  by  placing  the 
source  in  one  beam  of  the  interferometer.  In  the  present  paper,  we  explore  the 
efficiency  of  background  suppression  with  a  double-beam  i nterferometer  having 
coaligned  beams  and  tailored  MTFs.  This  new  technique  is  compared  to  the  back¬ 
ground  suppression  achieved  by  offsetting  the  two  beams. 

Consider  the  following  idealized  double-beam  i nterferometer  with  a  dielec¬ 
tric  beamsplitter  illustrated  in  Figure  1,  and  its  possible  outputs  illustrated 
in  Figure  2.  For  radiation  wi Lh  a  broad  spectral  distribution  B(a),  entering 
at  A  (beam  1),  the  output  of  detector  D  as  a  function  of  path  difference  x  is 

i  i 

given  by: 

a2 

D  (x) 

1  (beam  1) 

where  a  is  the  wave  number  and  o  and  a  are  the  spectral 

1  2 

limits  imposed  by  the  system  detector  or  a  filter. 


I  B  (a) 
Ja  1 


1  +  I 

2  2 


cos  2tfox 


da. 


(1) 


IGURE  1:  Schematic  of  double-beam  interferometer  with  a 
tailored  modulation  transfer  function  (MTF). 
Beams  of  radiation  entering  the  two  arms  of  the 
interferometer  pass  through  optical  components 
atAj  and  A  which  define  the  MTFs  for  the  two 
arms.  Mirrors  M  and  M  are  roof- type  retro- 
flectors.  The  dielectric  beamsplitter  is 
labelled  B.S. . 
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FIGURE  2:  Schematic  i nterferograms  for  a  dual-beam  interferometer,  recorded  by  a 
single  detector  (D^).  a)  radiation  of  integrated  intensity  B  entering 
arm  1  only;  b)  radiation  entering  arm  2  only;  c)  equal  fluxes  of  back¬ 
ground  radiation  entering  both  arms;  d)  equal  background  fluxes  and  a 
highly  unbalanced  source  flux  in  the  two  arms;  e)  an  equal  but  time- 
varying  background  flux  entering  both  arms;  and  f)  difference  in  inter- 
ferograms  from  detectors  0}  and  D2  with  equal  but  tine  varying  background 
fluxes  and  unbalanced  point  source  fluxes. 
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The  interferogram ,  D  ,  will  look  qualitatively  like  the  curve  in  Figure  2a; 
for  illustrative  purposes,  we  have  used  a  spectral  distribution  B(o)  which  is 
Gaussian  in  shape  and  centered  on  a  non-zero  frequency. 

Radiation  which  enters  at  A  (beam  2,  Fig.  1)  strikes  the  opposite  face 
of  the  main  interferometer  beamsplitter  (at  0).  Beam  2  is  assumed  to  be 
incoherent  with  respect  to  beam  1.  The  interferogram  from  beam  2,  as  recorded- 
by  detector  D  ,  is  then  complementary  to  the  i nterferogram  that  results  from  beam 
1  (see  Fig.  2b),  and  is  given  by: 


D(x)  =  f  B  (a)  \ 

1  (beam  2 )  Ja  1 


COS  2iraX  da 


(2) 


l 

The  interferogram  from  beam  2  takes  this  form  because  each  frequency  component 
undergoes  an  additional  n-  phase  shift  from  its  interactions  with  the  beam¬ 
splitter.  By  suitably  designing  the  instrument,  the  spectral  intensities  in 

r  °2 

beam  1  and  beam  2  can  be  made  equal.  Then  ^  ( x )  ( beams  1  and  2)  =  i  2^  da. 
If  B(a )  is  not  time  varying,  the  detector  output  at  D  is  a  constant.  From 


] 

the  symmetry  of  the  figure,  the  output  of  D  is  also  constant. 


2 

Consider  now  the  case  where  the  signals  introduced  into  the  two  arms 
of  the  interferometer  are  unbalanced.  Specifically,  let  the  signals  from  a 
point  source  in  the  field  of  view  be  unbalanced  (i.e.,  S  (o)  in  arm  1  and  S^(a) 
in  arm  2),  while  the  background  levels  are  still  quite  closely  balanced 
(B^o)  =  B  (u)  =  B(a)/2) .  Several  techniques  are  available  for  producing  the 
imbalance  in  source  signals  (while  leaving  the  background  levels  balanced),  and 
these  will  be  discussed  in  Section  2.2.  The  imbalance  in  target  signals  will  be 
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i. 


recorded  in  detector  Dj  as: 


D(x) 

!  (beams  1&2) 


B  +  \  (S,  +  S? )  +  \ 


(3) 


(OS  2-rroX  do. 


where  B  =  /  B(o)da  and  Sj  and  S?  are  similarly  defined.  This  intcrferogram  is 
shown  schematically  in  Figure  2d.  The  observed  modulation  amplitude  is  equiv¬ 
alent  to  signal  of  intensity  S[  -  S2  entering  only  one  arm  of  the  interferometer. 
This  modulating  signal  from  the  source  is  superposed  on  the  DC  level  provided 
by  the  source  plus  background  (see  Fig.  2d). 

The  preceding  assumes  the  background  is  not  fluctuating  with  time.  It  is 
also  possible,  however,  to  suppress  temporal  background  fluctuations.  Although 
such  fluctuations  would  occur  in  both  beams,  they  would  be  in  phase  and  result 
in  a  spurious  modulation  in  Dj  (Fig.  2e),  which  reduces  the  sensitivity  for  source 
detection.  We  note  that  this  limitation  can  be  overcome  by  making  use  of  the 
symmetry  of  the  interferometer  (Fig.  1): 


D(x) 

2  (beams  1&2) 


(4) 


B  +  \  (Si  +  S2)  + 


f 

[s 2  ( a )  -  S !  (a )  I 

cos  2tt ax  da 

L  J 

If  we  now  electronically  subtract  the  outputs  of  detectors  Dj  and  D? : 

D(x)  -  D(x)  =  /  FS^a)  -  S2(a)l  cos  2ttox  da.  The  contribution  from  equal 
l  2 

background  fluxes  in  the  two  arms  cancels,  while  the  source  modulations  are 
addi tive. 

2 . 2  Tailored  Modulation  Transfer  Functions 

The  background  suppression  in  the  i nterferogram  represented  by  equation 
3  is  based  on  the  assumption  that  one  can  obtain  a  significant  imbalance  in  the 
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source  signal  with  no  imbalance  in  the  background  signal.  As  a  simple  example, 
consider  a  point  source  viewed  against  a  uniform  background.  If  the  optical 
axes  of  beam  1  and  beam  2  are  slightly  offset,  then  the  two  images  of  the  beams 
will  be  displaced  in  the  image  plane.  For  an  appropriately  sized  detector 
no  suppression  of  the  point  source  i nterferograms  will  occur.  However,  dis¬ 
placement  does  not  affect  a  uniform  background  and  the  modulation  amplitude 
in  the  background  i nterferogram  will  be  suppressed.  The  problem  with  this  scheme 
is,  of  course,  that  the  background  may  not  be  uniform  from  point  to  point,  and 
its  suppression  will  therefore  be  incomplete. 

In  a  first  attempt  to  find  a  more  efficient  means  of  background  sup¬ 
pression,  we  considered  deliberate  introduction  of  a  small  amount  of  defocus 
into  one  of  two  coaligned  interferometer  arms.  Images  from  either  of  the  arms 
can  be  slightly  defocussed  on  the  detector  face  by  placing  a  weak  positive  or 
negative  element  at  A]  or  A2  (Fig.  1).  This  has  the  overall  effect  of  sup¬ 
pressing  the  modulation  amplitude  for  the  more  smoothly  varying  regions  of 
the  background,  while  not  suppressing  the  modulations  of  point-like  sources 
(see  Eq.  3).  Intuitively,  this  effect  can  be  understood  by  noting  that  a 
slightly  defocussed  image  of  a  smoothly  varying  background  does  not  differ 
much  from  the  correspondi ng  in-focus  image.  Thus,  when  the  in-focus  back¬ 
ground  image  from  one  beam  is  superimposed  in  coincidence  on  the  defocused 
image  from  the  other  beam,  the  i nterferogram  is  reduced  to  nearly  a  DC  level 
on  each  detector  element,  and  effective  background  suppression  occurs.  However, 
a  slightly  defocussed  image  of  a  point-like  source  differs  strongly  from  the 
corresponding  in-focus  image;  thus,  the  source  i nterferogram  is  highly  modulated 
when  the  two  images  are  superposed  on  the  detector  face. 
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Defousing  is  the  only  one  of  many  possible  methods  for  altering  the  response 
of  optical  systems  of  this  type;  the  two  MTFs  corresponding  to  the  two  arms  of 
the  interferometer  can  readily  be  tailored  to  satisfy  various  system  require¬ 
ments.  In  addition  to  defocusing,  some  of  the  other  techniques  which  can  be 
used  in  altering  MTF  pairs  are  listed  as  follows: 

1.  Introduction  of  controlled  amounts  of  an  aberration  or  some  special 
optical  component  into  the  system  (e.g.,  spherical  aberration,  a 
weak  axicon). 

2.  Adjustment  of  the  ratio  of  the  diameters  of  the  two  entrance  apertures. 

3.  Phase  apodization  of  an  entrance  aperture  by  the  introduction  of 
phase-shifting  annular  rings  or  disks. 

4.  Transmission  apodization  of  an  entrance  aperture  by  the  introduction 
of  annular  rings  or  disks  of  neutral  density  filters. 

5.  Controlled,  small  amplitude  mechanical  oscillation  of  an  optical 
component  in  the  system  at  a  frequency  high  compared  to  the  system 
response  time. 

The  above  techniques  can  be  applied  individually  or  in  combination  to  one  or  both 
arms  of  the  double-beam  interferometer  to  produce  changes  in  the  system  MTF. 

Consider  a  double-beam  interferometer  comprised  of  two  arms  each  having 
different  responses  to  various  spatial  frequencies.  When  the  two  beams  are 
combined  at  the  detector,  the  effect  is  that  of  a  spatial  filter  that  ideally 
responds  to  the  signal  from  the  point  source  being  studied,  while  suppressing 
the  background  signal  (see  Eq.  3).  The  effective  modulation  transfer  function 

(MTF)  of  the  double-beam  interferometer  is  the  difference  in  MTFs  of  the  two 
arms.  For  circularly  symmetric  beams,  the  effective  MTF  can  be  written  as 
(see  e.g. ,  ref.  9) : 


mrn 
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(5) 


MTF(K)  u  j*  I  I  ,(R)  -  I ; ( R )  I  J  ( 2 tt R K )  RdR 
o 

In  this  expression,  I  (R)  is  the  diffraction  pattern  as  a  function  of  radial 
distance  in  the  image  plane,  JQ  is  a  zero  order  Bessel  function,  and  K  is  a 
spatial  frequency.  We  emphasize  that  the  "subtraction"  process  discussed  here 
refers  only  to  the  modulation  amplitude  in  the  interferogram,  and  not  to  the 
DC  signals.  We  assume  here  that  the  wavelength  range  is  sufficiently  restricted 
(i.e.,  AX/ A  <  0.1)  that  the  MTF  may  be  assumed  to  be  a  constant  for  all  spectral 
wavelenqths  contributing  to  the  interferogram. 

For  purposes  of  illustration,  we  consider  here  tailored  pairs  of  MTFs 
that  result  when  one  of  the  following  four  types  of  aberrations  is  introduced 
into  the  system.  These  will  be  specified  by  the  aberration  function  6 ( p ) , 
where  p  is  the  distance  from  the  center  of  the  pupil  in  units  of  its  radius  a. 
The  aberrations  are:  1)  defocus  l<s(p)  =  Cp?J;  2)  n-phase  shift  L6  =  n  for 
p-  t  I;  3)  linear  aberration  (axicon)  i  cS ( p )  =  Cp'l;  and  4)  "wobbled  lens"  system. 
Figure  3  illustrates  the  tailoring  of  pairs  of  MTFs  in  the  ways  listed  above, 
with  the  goal  of  suppressing,  as  completely  as  possible,  the  system  response 
to  all  spatial  frequencies  below  about  one-fourth  the  cutoff  frequency 
(KfiWx  =  2a/x).  The  MTFs  were  numerically  computed  from  equation  5.  In 
Figures  3a,  b,  and  c  the  curve  labelled  "A"  represents  the  MTF  for  the  path 
of  the  interferometer  with  no  aberrations.  The  curves  labelled  "B"  are  the 
MTFs  for  the  specified  aberration.  The  heavy  curve,  labelled  "A-B",  is  the 
effective  MTF  of  the  system. 

At  low  spatial  frequencies,  the  MTFs  shown  in  Figures  3a,  b,  and  c  are 
characterized  by  MTF(K=0)  =  0  followed  by  an  approximately  linear  rise  with 
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FIGURE  3:  Modulation  transfer  functions  for  several  pairs  of  aberrated  lens 
systems.  In  a),  b) ,  and  c)  the  curve  labelled  A  is  the  MTF  for  an 
ideal  lens.  The  curve  labelled  B  is  the  MTF  corresponding  to  the 
aberration  which  has  been  introduced  into  one  of  the  arms.  The 
curve  A-B  represents  the  effective  MTF  of  the  interferometer.  The 
"wobbled  lens"  system  in  d)  is  described  in  the  text.  Kmax  is  the 
diffraction  limited  cut-off  frequency,  2a/A,  where  a  is  the  aper¬ 
ture  (beam)  radius  and  A  is  the  wavelength  of  the  radiation. 
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K.  Note  that  for  these  3  cases  there  is  no  spatial  frequency,  other  than 
K=0,  which  is  perfectly  filtered.  Of  the  MTFs  shown  in  Figure  3,  defocus  is 
the  most  convenient  for  use  in  the  double-beam  interferometer  because  the 
MTF  may  be  changed  without  the  need  for  fabrication  or  installation  of  new 
optical  components. 

As  an  example  of  a  more  efficient  spatial  filter,  we  have  devised  at 
least  one  system  with  MTF(K)  =  0  for  all  K  <  0.2K  .  This  system  consists  of 

a  lens  with  a  n-phase  shift  for  p  <  0.3,  and  another  lens  that  has  a  n-phase 
shift  for  p  0.98,  and  that  is  rotated  off-axis  (wobbled)  at  a  frequency 
rapid  compared  to  the  system  response  time.  The  lens  is  wobbled  in  such  a 
manner  as  to  cause  the  point  spread  function  in  the  image  plane  to  rotate  in 
a  circle  at  a  constant  angular  velocity.  The  angular  diameter  of  this  circular 
motion  was  set  at  3'/8a.  The  chosen  parameter  values  for  the  phase  shifts  and 
the  angular  extent  of  the  wobble  are  somewhat  arbitrary  in  that  other  combina¬ 
tions  can  also  yield  an  effective  spatial  filter.  The  computed  MTF  for  this 
system  is  shown  in  Figure  3d.  The  only  disadvantage  of  this  system  is  that  it 
is  somewhat  more  difficult  to  implemert  than  those  summarized  in  Figures  3a, 
b,  and  c. 

*► 

The  signal  S(R)  recorded  by  a  detector  in  the  focal  plane  array,  from 
a  source  located  at  a  position  R  relative  to  that  detector,  is  given  by  the 
convolution  of  the  system  point  spread  function,  the  source  shape  and  the 
detector  shape.  This  is  equivalent  to  the  inverse  Fourier  transform  of  the 
product  of  the  respective  transforms.: 

-> 

S(R)  “ J  1 opt(K)  's(K)  'd(K)  e1K'RdK  (6) 
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In  this  integral  t  ^  is  the  effective  MTF  of  the  optical  system,  ts  is  the 
transform  of  the  spatial  distribution  of  the  source  as  projected  on  the  detector 
and  is  the  transform  of  the  detector  shape. 

The  response  of  the  interferometer  to  the  background  radiation  will 
depend  critically  on  the  power  spectrum  of  the  background  spati.il  structure. 

The  rms  signal  from  the  background  (noise)  that  is  recorded  by  <i  given  focal 
plane  detector  element,  as  the  interferometer  is  pointed  in  various  directions, 
is  (cf .  Pried) 10 : 


(7) 


In  this  expression,  3 ( K)  is  the  two  dimensional  power  spectrum  of  the  spatial 
distribution  of  the  background. 

Ideally,  the  MTF  is  designed  to  optimize  the  signal  to  noise.  In  most 
situations  of  interest,  the  background  spectrum  has  strong  contributions  from 
low  spatial  frequencies  while  the  source  of  interest  is  close  to,  or  smaller 
than,  the  diffraction  limit  of  the  optical  system.  We  have  thus  attempted 
to  construct  MTFs  with  a  minimum  response  to  low  spatial  frequencies  while 
maximizing  the  high  frequency  response  (Fig.  3). 

We  have  evaluated  the  performance  (i.e.,  signal -to-noi se)  of  a  variety 
of  MTFs  for  several  general  classes  of  background  structure.  We  discuss,  as 
an  example,  backgrounds  with  spatial  structure  characterized  by  a  two  dimen¬ 
sional  power  spectrum  proportional  to  K"n.  It  can  easily  be  shown  that  a 
uniform  background  or  a  linearly  varying  background  is  completely  eliminated 
by  any  of  the  spatial  filters  considered. 


Consider  two  idealized  classes  of  MIT;  these  are  shown  as  an  inset  to 

Figure  4.  The  first  (Type  I)  has  i  ^  ( K ) =0  for  all  K  below  some  cutoff  spatial 

frequency  KQ.  In  the  second  filter  (Type  II),  ■;  increases  linearly  from 

K=0  (with  slope  a)  until  its  value  is  equal  to  that  of  an  ideal  lens,  from 

which  point  (also  denoted  K  )  it  matches  the  ideal  lens.  When  computing  the 

rms  noise  from  equation  (7),  one  generally  must  assume  a  low  frequency  cutoff 

(usually  set  by  the  size  of  the  field  of  view)  since  a  background  with  a  K"n 

spectrum  does  not  have  a  well  defined  variance.  In  our  calculations  we  have 

assumed  a  field  of  view  equal  to  1000  x  1000  resolution  elements;  the  results 

are  relatively  insensitive  to  this  choice.  For  simplicity  in  the  calculations 

we  have  also  a  ;sumed  circularly  shaped  pixels  with  an  angular  diameter  of 

'/2a.  The  signal -to-noi se  ratios  for  both  of  the  above  idealized  MTFs  have 

_2 

been  computed  for  a  range  of  background  power  spectra;  the  results  for  a  K 
spectrum  (see,  e.g.,  ref.  10  and  11)  are  shown  as  an  example  in  Figure  4. 

The  ordinate  is  the  ratio  of  signal -to-noise  of  the  specified  system  to  that 

for  an  ideal  lens.  The  optimum  cutoff  for  a  Type  I  filter  is  -0.7K .  The 

Type  II  filter  becomes  highly  effective  as  K  ->1  even  though  its  peak  response 
to  point  sources  becomes  very  small.  For  comparison,  the  signal -to-noise 
ratios  for  the  "wobbled  lens"  (Fig.  3d)  and  "defocus"  (Fig.  3a)  spatial  filters 
are  also  indicated  in  Figure  4. 

As  a  comparison  to  spatial  filters  which  utilize  coaligned  beams  and 

tailored  MTFs,  we  have  repeated  the  signal -to-noise  calculations  for  the 

_2 

case  of  an  interferometer  wth  two  offset  beams.  The  results  for  a  K  back¬ 
ground  power  spectrum  are  shown  in  Figure  4  as  a  function  of  the  offset  angle. 

The  background  suppression  is  optimized  for  an  offset  approximately  equal  to 
the  diffraction  limited  angular  resolution  of  the  instrument.  For  large  offsets. 
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FIGURE  4:  Sample  signal-to-noise  calculations  for  a  double-beam  interfer¬ 
ometer,  with  a  tailored  modulation  transfer  function.  The 
example  is  for  a  point-source  viewed  against  a  background  whose 
spatial  components  are  described  by  a  two-dimensional  power 
spectrum  K*2.  For  idealized  Type  I  and  Type  II  spatial  filters 
(see  text),  the  calculated  signal -to-noi se  is  plotted  vs.  the 
parameter  K0  (defined  in  the  inset).  The  limiting  case  of  Ko-*0 
represents  an  unfiltered  lens  system.  For  comparison,  the  signal- 
to-noise  calculations  for  a  double-beam  interferometer  with  offset 
beams  are  shown  as  a  function  of  the  off-set  angle  S  (upper 
abscissa;  6  is  expressed  in  units  of  60  =  A/2a).  The  signal-to- 
noise  for  the  special  case  of  a  "wobbled  lens"  system  (Fig.  3d) 
and  a  "defocused"  system  (Fig.  3a)  are  also  indicated  on  the 
figure  for  comparison. 
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the  background  suppression  becomes  ineffective  because  the  instrument  becomes 
sensitive  to  low  spatial  frequencies.  As  the  offset  approaches  zero,  the  signal- 
to-noise  vanishes  because  the  signal  from  the  source  is  itself  suppressed.  When 
optimized  for  a  K“^  background  power  spectrum,  this  type  of  spatial  filter 
is  about  a  factor  of  3  less  efficient  than  filters  designed  with  tailored  MTFs. 

We  have  evaluated  the  background  suppression  efficiency  for  several 
selected  filters  as  a  function  of  the  spectral  index  of  the  assumed  power- 
law  spectrum  (K  n)  which  describes  the  spatial  structure  in  the  background. 

In  all  cases,  the  spectrum  is  assumed  to  be  a  function  of  the  magnitude  of  K 
only.  The  results  are  shown  in  Figure  5  for  a  Type  I  filter  (K  =  0.20), 

Type  II  filter  (Kq  =  0.5),  "wobbled  lens"  system  (Fig.  3d),  defocused  system 
(fig.  3a),  and  an  interferometer  with  offset  beams  (offset  =  X/2a).  For  small 
spectral  indices  (n  <  1)  none  of  the  trial  spatial  filters  significantly 
enhances  the  signal-to-noise.  At  n=2  the  background  suppression  is  significant, 
with  the  "wobbled  lens"  system  yielding  a  factor  of  3  improvement  over  the 
offset  pointing  mode.  As  n  increases  (n  N  4)  there  are  substantial  differences 
in  performance  among  the  filters  of  up  to  two  orders  of  magnitude.  The  filters 
with  effectively  no  response  below  -0.2K  are  superior.  The  offset  pointing 
mode,  which  does  not  provide  a  good  "average"  of  the  background  around  the 
source,  is  the  least  effective. 


3. 


THE  VISIDYNE  DOUBLE -BEAM  INTERFEROMETER 


During  the  past  year,  Visidyne  has  designed  and  built  a  laboratory 
double-beam  i nterferometer- spectrometer  for  use  in  its  background  optical 
suppression  scheme.  The  working  i nterferometer  system  presently  features  a 
double-beam  interferometer  (2-cm  beam  aperture),  tailored  modulation  transfer 
function  for  the  two  beams,  servo-controlled  drive  for  the  moving  retroreflec- 
ting  mirror,  laser  retardation  reference  system,  PbS  detector  and  associated 
electronics,  16K  hardwired  fast  Fourier  transform  (FFT)  analyzer,  and  CRT 
plus  X-Y  plotter  for  displaying  interferograms  and  spectra  in  real  time.  The 
laboratory  facility  also  includes  a  test  setup  that  can  simulate  a  point 
source  immersed  in  a  controllable  sinusoidally  varying  background  at  infinity. 

The  general  layout  of  the  interferometer  laboratory  is  pictured  in 
Figure  6  and  shown  schematically  in  Figure  7.  The  interferometer,  the  data 
system  and  the  test  facility  are  described  in  detail  in  §  3.1,  3.2  and  3.3, 
respectively. 

3 . 1  The  Double-Beam  I nterferometer-Spectrometer 

The  Visidyne  laboratory  interferometer  is  a  double-beam  input  and 
dual-detector  output  design  with  CaF^  lenses  to  provide  a  tailored  modulation 
transfer  function  for  the  suppression  of  spatially  structured  background. 

The  interferometer  is  pictured  in  Figure  8  and  is  shown  schematically  in 
Figure  9a  (upper  level)  and  Figure  9b  (lower  level). 

3.1.1  Optical  Design 

Radiation  enters  the  system  through  either  a  flat  CaF^  compensator 
plate  (beam  1)  or  a  defocusing  CaF^  lens  (1  m  focal  length;  beam  2).  The 
beam  diameters  are  separately  variable  up  to  a  diameter  of  2  cm  using 


FIGURE  7:  Schematic  layout  of  i nterferometer  -  spectrometer  laboratory. 
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Side  view  of  interferometer. 


Top  view,  showing  beamsplitter,  retroflecting  mirrors,  detector  H, 
gas  cells,  and  laser  retardation  reference  system. 


FIGURE  8:  Double-beam  interferometer-spectrometer. 


FIGURE  9a:  Schematic  of  double-beam  interferometer-spectrometer  (upper  level). 
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FIGURE  9b:  Schematic  of  double-beam  interferometer-spectrometer  (lower  level). 
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adjustable  iris  aperture  stops.  Each  beam  is  then  incident,  from  opposing 
directions,  on  a  7.5  cm  diameter  beamc pi i tter- compensator  consisting  of  two 
parallel  plates  of  CaF^  between  which  is  sandwiched  a  thin  film  of  Si.  The 
transmission  of  the  beamsplitter  as  a  function  of  wavelength,  for  0°  inci¬ 
dence  angle,  is  shown  in  Figure  10.  The  transmitted  and  reflected  beams  are 
both  directed  to  roof- type  retrorefl ecting  mirrors.  One  of  the  two  retrore- 

flectors  is  movable;  both  displace  the  beam  down  to  the  lower  level  of  the 

o 

i nterferometer .  Roof- type  retrorefl ectors  were  selected  because  (1)  the 
beam  is  reflected  by  180°,  independent  of  small  tilts  of  the  retroreflectors 
about  two  of  its  three  axes  of  rotation,  and  (2)  the  use  of  roof  retroreflec¬ 
tors  allows  a  more  compact  i nterferometer  than  an  alternate  "cat's  eye"  de¬ 
sign.  The  mirrors  in  the  retrorefl ectors ,  as  well  as  all  other  reflecting 
surfaces  in  the  interferometer,  were  gold  coated.  The  i nterferometer  optics 
and  beamsplitter  have  been  designed  for  use  in  the  wavelength  range  1-10  urn 
In  the  lower  level  of  the  i nterferometer ,  the  beams  are  again 
partially  reflected  and  partially  transmitted  by  the  beamsplitter,  and  a  por¬ 
tion  of  the  radiation  that  enters  beam  1  arrives  at  detector  I.  The  comple¬ 
mentary  fraction  of  the  radiation  is  directed  to  detector  II.  The  detector 
consists  of  a  PbS  cell  onto  which  the  beam  is  focused  by  a  2.5  cm  diameter, 
f/2  CaF^  lens.  The  sensitive  area  of  the  detector  is  0.25  mm  x  0.25  mm. 
Ambient  temperature  PbS  detectors  were  used  for  these  measurements  because  of 
the  ease  with  which  detector  size,  a  BOSS  measurement  parameter,  could  be 
varied.  Each  detector  assembly  also  contains  a  spectral  filter  which  limits 
the  overall  bandpass  of  the  system.  The  spatial  background  suppression  data 
were  obtained  with  a  2.1  urn  to  2.5  urn  interference  filter.  This  bandwidth 
was  selected  so  as  to  minimize  atmospheric  absorption  effects  on  the 
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measurements.  Spectral  measurements  were  made  with  either  a  silicon  or 
germanium  window  in  place  of  the  filter. 

There  are  numerous  degrees  of  freedom  for  adjusting  and  aligning 
the  optical  components.  The  moving  retroreflector  is  mounted  on  a  micrometer 
adjustable  X-Y-Z  translation  stage.  The  two  mirrors  in  each  retrorefl ector 
can  be  separately  aligned.  The  orientations  of  each  of  the  other  mirrors  in 
the  interferometer  as  well  as  the  positions  of  the  beamsplitter,  the  fixed 
retroreflector  and  the  detectors  are  also  adjustable. 

As  the  retrorefl ecting  mirror  is  moved,  monochromatic  radiation 
that  enters  beam  1  and  arrives  at  detector  I  will  be  modulated  sinusoidally 
(see  §  2).  The  modulation  rate  is  proportional  to  the  frequency  of  the  in¬ 
cident  radiation  and  to  the  speed  of  the  mirror.  For  a  range  of  wavelengths 
entering  the  interferometer  in  beam  1,  a  maximum  in  inter  si ty  will  be  record¬ 
ed  in  detector  I  when  the  path  length  to  the  moving  retrc  ref  1  ector  equals  the 
path  length  to  the  fixed  mirror  (zero  path  difference).  (See  §  2  for  a 
detailed  discussion  of  the  optical  paths  through  the  interferometer.)  Radia¬ 
tion  which  enters  beam  2  and  arrives  at  detector  I  is  also  modulated  as  the 
retroreflector  moves.  If  the  two  beams  I  and  2  are  equal  (i.e.,  optical  axes 
viewing  the  same  direction,  balanced  intensities,  etc.)  the  modulations  will 
be  180°  out  of  phase  and  will  yield  a  DC  signal.  An  imbalance  in  the  two 
beams,  as  viewed  by  detector  I,  will  result  in  a  modulation  amplitude  propor¬ 
tional  to  the  imbalance. 

3.1.2  Mechanical  Design 

All  of  the  optical  components  of  the  interferometer  are  mounted  to 
an  aluminum  tooling  plate  of  dimensions  20"  x  36"  x  1"  (Figure  8).  The  tool¬ 
ing  plate  is  supported  above  a  granite  optical  bench  by  extruded  aluminum 
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beams.  Vibration  damping  of  the  table  is  provided  by  Firestone  air  mounts 
between  the  table  legs  and  the  laboratory  floor. 

The  retroreflector  drive  which  is  mounted  below  the  interferometer 
base  plate  is  shown  schematically  in  Figure  11.  The  drive,  based  on  a  design 
that  utilizes  a  precision  air-bearing  (Dover  Instruments,  Inc.),  a  precision 
lead  screw  (Precision  Thread  Grinding),  and  servocontrol led  motor  (Torque 
Systems),  has  proved  very  successful.  The  system  has  been  tested  for  a  mir¬ 
ror  travel  distance  of  1  cm  (see  §  4),  and  will  be  used  in  the  future  with  a 
10  cm  travel.  To  incorporate  the  10  cm  mirror  travel  requires  only  transla¬ 
tion  of  the  retrorefl ectors  to  other  mounting  holes  that  are  already  provided 
on  the  mounting  plate. 

The  speed  of  the  motor  is  controlled  by  a  servo- system  shown  schem¬ 
atically  in  Figure  12.  The  mirror  position  is  controlled  by  the  ramp  output 
from  a  waveform  generator.  The  errors  in  the  instantaneous  position  and 
speed  of  the  mirror  are  then  fed  back  into  the  servo-ampl i fier  to  maintain 
constant  mirror  velocity. 

Several  of  the  opto-mechanical  components  of  the  system  are  com¬ 
pletely  interchangeable  for  purposes  of  checking  the  balance  and  synmetry  of 
the  system.  These  include  the  detectors,  detector  mounts,  retrorefl ecti ng 
mirrors,  lens  mounts,  iris  diaphragms,  and  gas  cells. 

3.1.3  Laser  Reference 

A  He-Ne  laser  (\x6328  $)  is  mounted  to  the  interferometer  to  pro¬ 
vide  an  absolute  retardation  reference  (Figure  8).  The  laser  beam  traverses 
basically  the  same  optical  path  as  that  for  beam  2  with  the  exception  of  by¬ 
passing  the  first  mirror  (Ml,  Figure  9a),  and  returns  to  its  own  silicon 
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FIGURE  12:  Interferometer  mirror-drive  servo  schematic. 


detec tor-ampl ifier .  The  laser  beam  utilizes  one  corner  of  the  retrorefl  ect- 
ing  mirrors  and  thus  avoids  mixing  the  laser  light  wi  .h  the  signal  entering 
the  principal  detectors  of  the  i nterferometer .  Provisions  have  also  been 
made  to  monitor  the  position  of  the  opposite  corner  of  the  retrorefl ectors, 
during  a  scan,  to  ensure  that  there  are  no  spurious  angular  motions  of  the 
mirror  during  scans. 

3.1.4  Alignment  Procedure 

Initially,  the  beam  splitter  and  the  two  retrorefl ectors  are  align¬ 
ed  with  an  autocollimator.  The  system  is  next  brought  into  finer  alignment 
visually  using  both  monochromatic  and  white  light  sources  and  observing  the 
fringe  patterns.  Adjustments  are  made  until  the  fringe  patterns  of  concen¬ 
tric  rings  are  centered  on  the  optical  axis  of  the  interferometer. 

The  two  input  beams  are  then  co-aligned  by  adjusting  M^.  Fine 
alignment  adjustments  are  made  with  the  detector  in  place  and  both  beams  of 
the  interferometer  viewing  a  diffuse  light  source.  When  the  ratio  of  single¬ 
beam  to  double-beam  i nterferogram  amplitudes  at  zero  path  difference  is  maxi¬ 
mized,  the  i nterferometer  is  aligned.  In  practice  the  maximum  value  of  this 
ratio  was  approximately  100. 

3.1.5  I nterferometer  Configurations 

The  interferometer  can  be  used  in  three  different  modes,  (i)  The 
single-beam  mode  (with  one  of  the  two  beams  blocked)  is  used  to  obtain  spec¬ 
tra  of  remote  sources.  ( i i )  The  double-beam  mode  with  a  tailored  modulation 
transfer  function  is  used  to  obtain  spectra  of  remote  localized  sources  in 
the  presence  of  a  spatially  structured  background,  (iii)  The  double-beam  and 
single  beam  modes  can  be  used  to  measure  the  absorption  (and  transmission) 
spectra  of  laboratory  gas  samples.  In  this  configuration,  the  sample  is 
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placed  directly  in  one  of  the  beams  and  light  from  an  external  source  is 
viewed  through  the  sample.  The  gas  cells  shown  in  Figure  8  are  used  for  this 
mode . 

3.2  The  Data  System 

The  data  system  for  the  i nterferometer- spectrometer  is  shown  schem¬ 
atically  in  Figure  13;  some  of  the  key  elements  are  pictured  in  Figure  14.  As 
the  retrorefl ecti ng  mirror  is  driven  in  linear  motion,  the  signal  from  the 
detector  is  amplified  and  filtered  through  a  low-pass  antialiasing  filter 
having  a  120  db/octave  roll-off.  Simultaneously,  the  laser  reference  detec¬ 
tor  electronics  sense  positive-going  zero-crossings  and  trigger  a  CORDIC  16K 
FFT  analyzer  (Shumway  Bros;  design  rights  acquired  from  Computer  Information 
Systems,  Inc.).  The  FFT  analyzer  then  samples  the  filtered  i nterferometer 
detector  output,  A/0  converts  the  signal  to  a  12  bit  word  and  stores  the  re¬ 
sult  as  one  input  point  of  a  16K  point  i nt.er ferogram.  The  analyzer  has  the 
capability  of  transforming  32K  words  per  second. 

While  the  i nterferogram  is  being  read  into  the  CORDIC  16K  FFT  anal¬ 
yzer,  the  data  are  transformed  and  an  8  K-point  amplitude  spectrum  is  output¬ 
ted  for  display  or  plotting.  Real  or  imaginary  transforms  can  both  be 
displayed.  The  analyzer  has  the  capability  of  performing  a  running  average 
of  output  spectra  for  a  range  of  selectable  averaging  times.  At  the  present 
time,  a  full  double-sided  i nterferogram  is  inputted;  in  the  future,  a  white 
light  reference  source  will  be  incorporated  into  the  i nterferometer  to  pro¬ 
vide  a  zero  path  difference  reference  that  will  enable  the  analyzer  to  pro¬ 
cess  single-sided  i nterferograms .  Because  the  phases  of  the  various 
components  in  the  i nterferogram  are  a  very  slowly  varying  function  of 
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FIGURE  13:  Double-beam  interferometer  laboratory  data  system. 


frequency,  one  can  use  this  information  to  recover  nearly  N  spectral  points 
from  an  N  point  i nterferogram' '  . 

Spectra  and  i nterferograms  are  displayed  in  several  ways.  Most  of 
the  spectra  shown  in  s  4  were  produced  directly  on  the  X  -  Y  plotter  (Figs. 

13  and  14);  others  were  computer-plotted  using  the  analyzer  digital  data 
output.  Inter ferograms  are  usually  displayed  on  a  CRT  screen  (see  Fig.  17). 

The  lb  K  CORDIC  analyzer  can  also  be  used  to  perform  much  longer 
transforms  such  as  would  be  obtained  from  data  with  0.1  cm  *  resolution.  A 
long  data  train  (e.g.,  2.6  x  10'  points)  would  be  read  into  the  computer  and 
stored.  The  data  train  is  multiplied  by  a  sinusoid  of  a  particular  frequency 
v0  and  then  summed  into  coarser  bin  sizes  (e.g.,  lb  of  the  original  i nterfero¬ 
gram  bins).  The  compressed  data  train  is  then  transformed  with  the  CORDIC 
analyzer  in  the  usual  way.  The  result  is  a  li>K  segment  of  the  high  resolu¬ 
tion  spectrum  starting  at  vQ.  This  procedure  is  repeated  until  the  entire 
transform  is  complete.  Ihis  method  is  discussed  in  detail  in  reference  13. 

3.3  The  Test  System 

A  special-purpose  test  and  calibration  system,  to  measure  the  spec¬ 
tral  and  spatial  filtering  properties  of  the  interferometer,  has  been 
designed  and  constructed.  The  system  is  pictured  in  Figure  15  and  shown 
schematically  in  1  igure  16.  A  "point  source"  is  produced  by  a  pinhole,  illum¬ 
inated  by  a  tungsten-halogen  lamp,  in  the  focal  plane  of  an  8"  diameter  f/5 
collimating  mirror.  The  mirror  views  the  pinhole  by  reflection  from  an  ellip¬ 
tical  flat  (to  remove  the  light  sources  from  the  optical  path  to  the  interfer¬ 
ometer),  and  a  beamsplitter  that  allows  a  spatially  structured  source  to  be 
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Side  view  showing  beamsplitter,  chopper  wheel  assembly  and  focal 
plane  light  sources. 


Front  view  showing  f/5  primary  mirror,  elliptical  flat  and  focal 
plane  light  sources. 


FIGURE  15:  Laboratory  test  system 
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FIGURE  16:  Schematic  of  laboratory  test  system. 
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superimposed  in  the  focal  plane  (Figure  16).  Spectral  structure  can  be  added 
to  either  the  point  source  or  the  spatially  structured  source  with  filters  of 
a  variety  of  materials. 

The  second  source  in  the  focal  plane  of  the  mirror  is  designed  to 
measure  the  response  of  the  i nterferometer  to  various  spatial  frequencies  in 
a  spatially  structured  background.  At  present,  this  source  consists  of  a 
pair  of  Ronchi  rulings  whose  Moir£  pattern  can  be  described  as  a  two  dimen¬ 
sional  plane  wave  with  a  triangular  waveform.  The  spatial  frequency  of  the 
wave  can  be  varied  continuously  from  DC  to  high  frequencies  above  the  re¬ 
sponse  of  the  instrument.  The  phase  of  the  wave  can  also  be  varied  either  by 
adjusting  the  relative  positions  of  the  two  rulings,  or  by  a  translation 
stage  that  shifts  the  pair  of  rulings.  The  Moir<?  pattern  is  partially  trans¬ 
mitted  through  a  beamsplitter  where  it  is  subsequently  collimated  by  the  f/5 
mirror  and  projected  in  the  direction  of  the  i nterferometer . 

Both  beams  of  the  interferometer  thus  view  identical  "scenes"  con¬ 
sisting  of  a  point  source  superimposed  on  a  plane  wave  of  variable  spatial 
frequency;  these  two  components  may  each  have  different  spectra.  The  results 
of  tests  with  this  apparatus  are  presented  in  §  4. 
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4. 


LABORATORY  TESTS  AND  DEMONSTRATIONS 


The  Vi  si  dyne  double-beam  i nterferometer- spectrometer  was  designed 
especially  to  enhance  the  detectability  of  small  localized  sources  viewed 
against  a  spatially  structured  background  in  the  infrared  region  of  the  spec¬ 
trum.  The  spectral  information  inherent  in  the  interferometer  also  allows 
for  additional  enhancement  of  the  signal-to-noise  if  the  localized  source  has 
a  different  spectral  signature  than  that  of  the  background.  The  technique  of 
using  a  double-beam  interferometer,  with  a  tailored  modulation  transfer  func¬ 
tion  to  suppress  a  spatially  extended  background,  while  responding  well  to 
localized  sources,  is  discussed  in  detail  in  §2.  The  design  and  fabrication 
of  the  laboratory  double-beam  interferometer-spectrometer  and  all  its  subcom¬ 
ponents  are  discussed  in  §  3.  Also  described  in  §  3  is  the  laboratory  test 
apparatus  for  measuring  the  response  of  the  interferometer  to  various  spatial 
and  spectral  frequencies. 

In  this  section  we  describe  measurements  of  the  spatial  filtering 
properties  of  the  interferometer  {§  4.1)  In  these  tests,  we  consider  only 
the  total  response,  i.e.,  the  integral  of  the  observed  spectrum  or,  approxi¬ 
mately,  the  amplitude  of  the  zero  path  difference  fringe  in  the  interfero- 
gram.  In  §  4.2  we  examine  the  background  suppression  as  a  function  of  both 
spatial  and  spectral  frequency.  Finally,  in  §  4.3  some  higher  resolution 
spectra  (at  1  cm"^  resolution)  are  presented  as  a  demonstration  that  an  inter¬ 
ferometer  operating  as  a  spatial  filter  need  not  sacrifice  frequency  resolu¬ 
tion.  (The  drive  in  our  interferometer  (see  §  3)  is,  however,  capable  of 
yielding  0.1  cm  ^  resolution,  and  will  be  so  tested  in  the  near  future.) 

The  prinicipal  raw  data  from  a  given  scan  of  the  interferometer  is 
an  i nterferogram.  Sample  i nterferograms  are  illustrated  in  Figure  17.  The 
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interferogram  from  beam  1  (focused  beam;  see  Fig.  1)  viewing  a  diffuse  back¬ 
ground  source  is  shown  in  Figure  17a,  and  th<>  correspondi ng  interferogram 
from  beam  2  (defocused  beam)  is  shown  in  17b.  When  both  beams  are  allowed  to 
combine  at  the  detector,  the  interferogram  is  as  shown  in  17c.  Because  the 
background  is  diffuse,  both  beams  have  a  nearly  identical  response;  however, 
the  interferograms  are  out  of  phase  with  respect  to  each  other  by  180°  (see  § 
2,  Fig.  2).  The  signal  is  thus  reduced  to  nearly  a  DC  level,  and  the  back¬ 
ground  is  thereby  suppressed. 

When  a  point  source  is  viewed  against  a  diffuse  background,  the 
results  are  quite  different.  The  interferogram  from  beam  1  is  shown  in  Fig¬ 
ure  17 d ,  and  that  from  beam  2  is  shown  in  Figure  17e.  In  the  douole-beam 
mode  (Fig.  1 7 f ) ,  the  diffuse  background  is  suppressed  as  discussed  above;  the 
interferogram  from  the  relatively  weaker  point  source  is  quite  discernible 
(note  the  differently  shaped  envelope).  This  is  a  result  of  the  fact  that 
beam  1  and  beam  2  respond  very  differently  to  a  point  source,  and  the  inter¬ 
ferogram  is  therefore  not  suppressed  (i.e.,  not  reduced  to  a  DC  level). 

In  the  remainder  of  this  section,  we  show  only  the  spectra,  which 
are  computed  as  Fourier  transforms  of  the  interferograms .  The  raw  data,  how¬ 
ever,  in  each  case  are  interferograms  of  the  kind  shown  in  Figure  17. 

4.1  Spatial  Filtering  Properties 

The  point  spread  function  of  beam  1  (focused  beam)  and  the  point 
spread  function  of  beam  2  (defocused  beam)  were  measured  separately.  The 
point  light  source  in  the  laboratory  test  apparatus  (Fig.  16)  was  used  with  a 
chopper  wheel  to  modulate  the  signal.  (The  retroreflectors  were  at  rest 
during  the  test,  and  hence  there  was  no  modulation  of  the  signal  due  to  vary- 
i ng  path  differences.)  The  measured  point  spread  functions  for  the  two  beams 
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are  shown  in  Figure  18;  these  include  the  effects  of  the  CaF^  defocusing 
lens,  the  detector  lens  and  its  aberrations,  and  the  finite  size  of  the  de¬ 
tector.  The  ratio  of  peak  responses  for  the  focused  to  defocused  beams  is 
~  3.3.  We  emphasize,  however,  that  the  "areas"  under  both  curves  [i.e., 
j  I  (  G  )2  tt  odo  where  1  (  o)  is  the  appropriate  point  spread  function]  are  equal. 

The  modulation  transfer  function  (MTF)  for  each  of  the  two  beams 
(1  and  2)  as  well  as  the  effective  MTF  for  the  system  (1  minus  2)  were  com¬ 
puted  by  taking  a  Bessel  transform  of  the  point  spread  functions  shown  in 
Figure  18.  Specifically,  we  carried  out  the  following  integral: 

MTF(K)  «  f  1(e)  J(2„Ke)  2uede  (8) 

where  J  is  a  zero  order  Bessel  function  and  K  is  the  spatial  frequency.  The 
results  are  shown  in  Figure  19.  Note  that,  as  expected,  the  MTF  for  the  de¬ 
focused  beam  falls  off  at  the  higher  spatial  frequencies  considerably  faster 
than  the  MTF  for  the  focused  beam.  The  system  MTF,  (applicable  when  the  in¬ 
terferometer  is  used  in  the  double-beam  mode)  has  a  value  of  0  at  K=0, 
exhibits  a  nearly  linear  rise  for  K  ~  0.5  cyc/deg,  peaks  at  ~  0.7  cyc/deg, 
and  is  followed  by  a  decay  to  zero  at  K  ~  2.5  cyc/deg.  Thus,  the  present 
optics  yield  the  greatest  sensitivity  for  sources  with  an  angular  extent  of 
„  0.4°-1.2°. 

The  sensitivity  of  the  double- beam  interferometer  to  various  spa¬ 
tial  frequencies  was  directly  measured  with  the  laboratory  test  source  discus¬ 
sed  in  §  3  (Fig.  16).  The  pair  of  Ronchi  rulings  in  this  source  was  adjusted 
to  yield  a  Moire  pattern  with  a  specific  spatial  frequency  (as  well  as  the 
higher  harmonics  associated  with  the  triangular  wave  form).  The  phase  of  the 
Moire  pattern  was  adjusted  relative  to  the  optical  axis  of  the  interferometer 
until  a  maximum  in  intensity  in  the  central  fringe  of  the 
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FIGURE  18:  Point  spread  functions  for  the  beams  with  the  CaF2  defocusing  lens 

and  the  CaF2  compensator  plate.  In  each  case,  the  spatial  resolution 
of  the  detector  is  included. 
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interferogram  was  recorded.  At  this  setting,  the  interferogram  and  the  trans¬ 
formed  spectrum  were  recorded  for  both  the  single-beam  and  double-beam  modes. 
The  phase  of  the  wave  was  then  shifted  until  a  minimum  in  intensity  was 
obtained,  and  data  were  similarly  recorded.  This  procedure  was  repeated  for 
a  number  of  different  values  oF  the  spatial  frequency  K.  The  response  of  the 
interferometer  (defined  as  the  spectrum  integrated  over  all  frequencies)  in 
the  double-beam  mode  was  then  normalized  to  the  response  of  the  focused  beam 
alone.  The  results  as  a  function  of  K  are  plotted  in  Figure  20;  the  data 
points  rise  to  a  pronounced  maximum  at  K  T  0.9  cyc/deg. 

In  the  above  test,  the  measured  response  of  the  i nterferometer  to 
different  spatial  frequencies  is  closely  related  to  the  MTF  of  the  system. 

To  see  this,  we  note  that  the  theoretical  response  of  the  i nterferometer  to  a 
two-dimensional  source  with  spatial  distribution  S(x,y)  and  transform 
S(kx,k.y)  is  (ref.  10): 

D(x,y)«  ^IS(kx,ky)  MTF(kx,ky)  elk*x+1V  dkxdky,  (9) 

(cf.  equation  6,  §  2).  For  the  triangular  waveform  produced  by  the  Ronchi 
rul i ngs, 

S(x,y)  =  SQ[1  +  8/*2  £  n"2  cos(nkQx)]/2  and  (10) 

S(kx,k  )  =  Sq1t  U(kx)  +  4 /it  £  n'2{6(kx-nkQ)+  6(kx+nkQ)}]  «( ky ) 

odd 

where  kQ  is  the  fundamental  spatial  frequency  of  the  Moir£  pattern.  There¬ 
fore  D(x,y)  can  be  computed  analytically: 

D( x,y)  =  SQ  [MTF  (0)  +  8/tt2£  n'2MTF(nkQ)  cos(nkQx)]  .  (11) 

odd 

The  ratio  of  the  maximum  response  in  the  double  beam  mode  to  the  single  beam 
mode  is  then: 

RU0)  =  MTFA_B(0)+8A2Zn-2MTFAB(nko) 

MTFa(0)  +8/tt2  X  n"2MTFA(nkQ)  (1?) 


51 


DOUBLE  BEAM /SINGLE  BEAM  RESPONSE 
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FIGURE  20:  Ratio  of  the  double-beam  response  of  the  interferometer  to  the  single¬ 
beam  response,  as  a  function  of  the  spatial  frequency  of  the  Moirt; 
pattern.  The  solid  and  dashed  curves  are  the  expected  ratios  (eq.  12) 
for  an  ideal  i nterferometer  with  95%  and  100%  suppression  of  the  DC 
background,  respectively. 
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where  MTFA  and  MTF^  B  are  the  MTFs  for  beam  1  and  beams  1  minus  2,  respect¬ 
ively.  For  the  case  where  MTFA (0)=0  (perfectly  balanced  beams),  and  for 
those  values  of  nkQ  where  MTFA( nkQ) <<  MTFA(0): 

3(k0):  8/,2X  n'2MTFA  B(nk0)  :  8/n2  MTFA_B(kQ)  .  (13) 

Odd 

Thus,  by  means  of  the  above  computations,  the  measured  ratio,  R,  of 
double-beam  to  single-beam  response  as  a  function  of  spatial  frequency  can  be 
directly  interpreted  in  terms  of  the  MTF  of  the  interferometer.  The  dashed 
curve  in  Figure  20  shows  the  expected  ratio  R(k)  computed  from  equation  (12) 
for  the  case  where  the  DC  component  in  the  spatially  structured  background  is 
perfectly  suppressed  (i.e.,  perfectly  balanced  beams  with  MTF ( 0 )  =0  ).  The 
solid  curve  is  computed  from  equation  (12)  for  an  imbalance  of  5,  between  the 
two  beams.  Note  that  the  data  points  qualitatively  follow  the  latter  curve. 
The  qualitative  differences  could  be  due  to  a  combination  of  the  following 
causes:  (i)  imperfect  alignment  between  the  two  beams,  (ii)  non-uniform 
illumination  of  the  Ronchi  rulings,  and  (iii)  edge  effects  due  to  the  finite 
size  of  the  test  pattern. 

In  the  future,  we  plan  to  replace  the  simple  fringe-pattern  source 
in  the  test  apparatus  by  more  complicated  spatial  structures  (e.g.,  aerial 
photographs)  that  contain  an  entire  spectrum  of  spatial  frequencies.  For 
diagnostic  purposes,  however,  the  present  use  of  Ronchi  rulings  permits  a 
better  understanding  of  the  instrument  response  to  specific  frequencies. 

4 . 2  Combined  S patial an d  Spectral  Filtering  Propert ies 

In  the  previous  section  (§  4.1)  we  described  the  spatial  filtering 
properties  of  the  i nter ferometer ,  and  in  particular  its  capability  for  sup¬ 
pressing  low  spatial  frequencies.  In  this  section  we  describe  how  the  back 
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ground  suppression  of  the  i nterferomter-spectrometer  is  enhanced  by  utilizing 
the  spectral  information  inherent  in  such  a  device. 

A  Moir£  pattern  of  a  given  spatial  frequency  was  produced  with  the 
Ronchi  rulings,  and  a  relatively  low-intensity  point  source  was  superimposed 
by  illuminating  the  pinhole  on  the  complementary  side  of  the  beamsplitter.  A 
spectral  filter  composed  of  1.5  rain  thick  polycarbonate  was  placed  in  front  of 
the  point  source.  Thus,  the  "background"  had  both  a  different  spectral  as 
well  as  spatial  composition  than  the  weaker  point  source  which  was  to  be 
detected. 

For  each  of  several  different  spatial  frequencies  the  i nterferogram 
and  the  transformed  spectrum  were  recorded  in  both  the  single-beam  and 
double-beam  modes.  The  spectra  obtained  for  three  sample  spatial  frequencies 
(  0,  1.1,  and  7  cyc/deg)  are  shown  in  Figures  21,  22,  and  23. 

The  response  to  a  uniform  background  (K=0)  with  point  source  super¬ 
posed  is  depicted  in  Figure  21.  Spectra  for  the  point  source  plus  background 
are  shown  for  the  double-beam  and  single  beam  modes;  a  double-beam  spectrum 
of  the  background  alone  is  also  shown  for  comparison.  The  overall  spectral 
shape  (~  4000-5000  cm-*)  exhibited  in  the  single-beam  mode  results  from  the 
transmission  properties  of  an  interference  filter  located  in  front  of  the 
detector.  In  tie  double-beam  mode,  nearly  all  of  the  background  radiation  is 
suppressed,  while  the  spectrum  of  the  polycarbonate  filter  placed  in  front  of 
the  point  sourca  is  clearly  revealed.  By  contrast,  about  75A  of  the  inten¬ 
sity  of  the  point  source  is  retained  in  the  double-beam  mode  because  the 
source  has  contributions  from  all  spatial  frequencies.  (The  baseline  of  this 
spectrum  has  bean  shifted  upward  in  order  to  separate  the  curves.)  The 
double-beam  spectrum  without  the  point  source  shows  how  completely  the  back¬ 
ground  is  supprassed.  For  all  but  a  small  range  of  frequencies  (~  4100-4250 
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SOURCE  PLUS  BACKGROUND  SPECTRUM 


FIGURE  22:  Single  and  double-beam  spectra  of  a  point  source  superposed  on  a 
background  with  spatial  frequency  1.1  cycles  per  degree.  Shown 
for  comparison  is  a  double-beam  spectrum  of  the  background  alone. 
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FIGURE  23:  Single  and  double-beam  spectra  of  a  point  source  superposed  on 

a  background  with  a  high  spatial  frequency  of  7  cycles  per  degree. 
Shown  for  comparison  is  a  double-beam  spectrum  of  the  background 
alone  to  demonstrate  the  high  degree  of  suppression. 
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cm"1)  the  suppression  ratio  is  ~  100:1.  For  a  point  source  whose  spectrum  is 
characterized  by  a  narrow  line,  or  a  set  of  narrow  lines,  the  enhancement  of 
signal-to-background  over  a  simple  photometer  could  easily  exceed  10^:1. 

The  background  was  next  changed  to  a  fringe  pattern  of  frequency 
~  1.1  cyc/deg.  This  spatial  frequency  is  in  the  vicinity  of  the  peak  of  the 
MTF  curve  for  the  double-beam  mode  (Fig.  19).  The  background  at  this  fre¬ 
quency  is  therefore  i nccmpl etely  suppressed.  The  point  source  spectrum,  how¬ 
ever,  is  still  quite  pronounced  above  the  background  in  the  double-beam  mode. 

At  high  spatial  frequencies  (~  3  cyc/deg  with  the  present  MTF),  the 
interferometer- spectrometer  responds  in  virtually  the  same  way  as  for  a  DC 
background  (Fig.  19).  A  comparison  of  Figures  23  and  21  verifies  this  concl¬ 
usion.  The  basic  difference  between  the  two  figures  is  a  change  in  overall 
intensity  of  the  background  signal  which  results  from  the  different  DC  levels 
associated  with  the  Ronchi  rulings  in  different  orientations. 

4.3  Spectral  Tests 

There  are  several  ways  in  which  the  spectral  information,  inherent 
to  a  double-beam  i nter ferometer- spectrometer ,  can  be  utilized  to  enhance  the 
visibility  of  a  source  by  suppressing  unwanted  backgrounds.  First,  one  beam 
of  the  interferometer  is  aimed  at  the  source  of  interest  while  the  other  is 
offset  so  as  to  exclude  the  source.  In  this  way,  spatially  diffuse  back¬ 
ground  components  are  suppressed  regardless  of  their  spectra,  while  the 
higher  spatial  frequency  components  in  the  background  can  be  distinguished  by 
their  spectra.  Second,  different  spatial  frequency  responses  (MTF)  can  be 
incorporated  into  the  two  beams,  both  of  which  have  the  same  view  direction. 
Again,  spatially  diffuse  background  components  are  suppressed  regardless  of 
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their  spectra,  while  the  higher  spatial  frequency  components  in  the  back¬ 
ground  can  be  distinguished  by  their  spectra,  (see  §  4.2)  Third,  the  source 
may  be  placed  directly  into  one  of  the  two  beams  of  the  i nterferometer .  All 
spectral  components  common  to  both  beams  (i.e.,  most  non- source  background) 
will  be  suppressed.  The  latter  technique  is  not  directly  appplicable  to 
remote  sensing  problems,  but  is  the  most  convenient  laboratory  demonstration 
of  the  spectral  capabilities  of  the  interferometer. 

For  these  tests,  we  made  use  of  the  interferometer  configuration 
shown  in  Figure  24.  A  cell,  containing  a  gas  whose  absorption  spectrum  is  to 
be  measured,  is  placed  directly  in  the  path  of  one  of  the  beams.  An  identi¬ 
cal  cell  (not  shown  in  Figure  24)  is  placed  in  the  other  beam  to  null 
any  effects  due  to  the  first  cell.  The  single-beam  and  double-beam  spectra 
for  the  case  of  air  in  both  cells  are  shown  in  Figure  25.  The  overall  spectral 
shape  for  the  single-beam  response  is  due  to  the  silicon  filter  (in  front 
of  the  detector)  at  short  wavelengths,  and  the  PbS  detector  at  long  wave¬ 
lengths;  some  of  the  absorption  features  are  due  to  water  vapor  in,  the  air. 
Note  the  excellent  cancellation  of  the  signal  in  the  double-beam  mode,  as 
expected  (see  §  2,  Figure  2c).  The  spectral  resolution  here  is  10  cm_l. 

Pure  NH^  vapor  at  1  atmosphere  pressure  is  now  introduced  into  one 
of  the  cells.  The  spectrum  taken  in  the  single-beam  mode  is  shown  in  Figure 
26.  This  spectrum  represents  approximately  the  transmission  of  NH^  vapor 
multiplied  by  the  smoothly  varying  sensitivity  function  of  the  interferom¬ 
eter. 

Most  effects  of  the  filter,  the  water  vapor  in  the  air,  and  any 
effects  other  than  absorption  by  NH^  are  easily  and  effectively  removed  by  a 
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spectrum  produced  in  the  double-beam  mode.  To  see  this,  let  the  ideal  spec¬ 
trum  recorded  in  the  single-beam  mode  oe  S(0)  (this  includes  the  source  spec¬ 
trum  and  the  spectral  response  of  the  interferometer),  and  define  Tg( a)  and 
Tc ( a )  as  the  transmission  of  the  test  gas,  and  the  transmission  of  all  other 
matter  (e.g.,  water  vapor,  gas  cell,  etc.)  between  the  light  source  and  the 
detector,  respectively.  Then  the  response  in  the  single-beam  mode  is: 

Rsb(a)=S(a)Tc(c)Tg(a)  (14) 

while  the  response  in  the  double-beam  mode  is: 

Rdb(a)=S(o)Tc(o)-S(o)Tc(o)Tg(o)  (15) 

=S(a)Tc(a)[l-Tg(a)] 

=S(a)Tc(a)Ag(o) 

where  Ag(o)  is  the  absorption  of  the  test  gas.  Thus,  the  double-beam  re¬ 
sponse  is  the  absorption  of  the  gas  multiplied  by  the  source  spectrum  and  the 
absorption  spectra  of  other  contaminants.  The  advantage  of  this  presentation 
is  that  all  spectral  features  that  are  unrelated  to  the  gas,  and  that  are 
away  from  regions  of  significant  absorption  by  the  gas,  are  suppressed. 

The  results  for  NH^  absorption  (double-beam  mode)  are  shown  in 
Figure  27.  The  prominent  ( +  v-j)  and  ( \»^  +  v^)  vibration/rotation  bands 
(ref.  14)  at  '4450  cm  1  and  '5050  cm  respectively,  with  fine  structure 
down  to  10  cm’1  are  clearly  evident.  The  absorption  and  transmission  spectra 
are  superposed  in  Figure  28  to  emphasize  the  one-to-one  correspondence  of  the 
features  and  to  demonstrate  their  reality. 

A  similar  experiment  was  performed  with  a  1.5  mm  thick  polycarbon¬ 
ate  absorber  instead  of  the  gas  cell.  The  superposed  transmission  and  absorp¬ 
tion  spectra  are  shown  in  Figure  29.  A  broad  absorption  feature  for 
wavelengths  greater  than  '2.2  um  (o  ~  4500  cm’1)  is  clearly  evident. 
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FIGURE  27:  Double-beam  absorption  spectrum  of  NFL  vapor  in  the  region  of 
the  (\>2  +  v^)  and  (v^  +  v.)  absorption  bands. 
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FIGURE  28:  Double-beam  and  single-beam  spectra  of  NH3  vapor,  superimposed 
for  comparison. 
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FIGURE  29: 


Double-beam  and  single-beam  spectra  of  a  1.5  mn  thick 
polycarbonate  filter. 


Another  experiment,  with  NF^  in  one  of  the  gas  cells,  was  perform¬ 
ed  with  a  spectral  resolution  of  1  cm  1 ,  The  absorption  spectrum,  taken  in 
the  double-beam  mode,  of  the  (  Vp  +  band  centered  near  4450  cin’1  is  shown 
in  Figure  30.  Fine  rotation  structure,  down  to  1  cm’1,  within  the  band  is 
clearly  evident.  An  NH^  absorption  spectrum  (double-beam)  near  the  same  4450 
cm  1  band  is  superposed  on  the  compl ementary  transmission  spectrum  (single¬ 
beam)  in  Figure  31.  The  one-to-one  correspondence  of  the  fine  spectral 
detail  again  demonstrates  the  reality  of  these  features. 

The  ability  of  the  instrument  to  suppress  spectral  structure  in 
the  background  (away  from  the  absorption  features  of  the  test  gas;  e.q.  15) 
is  demonstrated  in  Figure  32.  The  upper  trace  is  the  single-beam  transmis¬ 
sion  spectrum  of  NH^  in  the  vicinty  of  the  ( +  v^)  band.  Also,  clearly 
present  is  the  ( ^  +  V3)  water  vapor  feature  (ref.  14)  centered  near  5330 
cm’1.  The  lower  trace  is  the  absorption  spectrum  of  NH^  ( doubl e-beam  mode) . 
Note  the  nearly  exact  correspondence,  in  both  traces,  of  the  NH^  features, 
while  the  water  vapor  features  above  5200  cm’1  are  completely  absent  in  the 
double-beam  spectrum.  The  ambient  water  vapor  is  common  to  both  beams  and  is 
thereby  eliminated.  In  remote  sensing  applications,  this  is  a  powerful  tool 
for  removing  the  spectral  contributions  from  any  background  sources  that  lie 
along  the  1 ine-of-sight  to  the  object  being  studied. 

Finally,  for  completeness,  we  show  the  1  cm’1  resolution  transmis¬ 
sion  spectra  for  laboratory  air  containing  water  vapor  in  the  and  bands 
centered  near  3700  cm_i  (Figure  33)  and  the  ( )  band  centered  at  5330 
cm’1  (Figure  34). 
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FIGURE  32:  Superposition  of  double-beam  (lower  trace)  and  single-beam  (upper  trace) 
spectra  of  NH3  vapor  in  the  vicinity  of  the  (.3  +  .4)  band.  Note  the 
unsuppressed  spectrum  of  H2O  vapor  (  -5200  cnrl]  in  the  single-beam  spectrum. 


X  ( microns) 

.80  2.75  2.70  2.65  2.60  2 

TTTTTTTTT  TT  n  I  I  I  I  I  TTT  TT  f  1 

WATER  VAPOR 


3600  3700  3800  390 

<r  ( cm'1) 


FIGURE  33:  Single  beam  transmission  spectrum  of  the  H„0  vapor  in  the 
vicinity  of  the  v,  and  v,  band?  (1  cm-1  resolution). 
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5.  CHARACTERISTICS  AND  APPLICATIONS  OF  DOUBLE-BEAM  INTERFEROMETRY 


In  this  report  we  have  demonstrated  both  empirically  and  computationally 
that  the  double-beam  interferometer-spectrometer  is  a  highly  efficient  tool 
for  the  remote  sensing  of  localized  sources  in  the  presence  of  strong  spa¬ 
tially  structured  background  in  the  IR  region  of  the  spectrum  (S  2  and  4). 

The  Visidyne  double-beam  interferometer  with  tailored  modulation  transfer 
functions  (§  2  )  utilizes  both  spatial  and  spectral  filtering  ( §  4)  to 
achieve  enchanced  sensitivity  to  localized  sources.  In  the  double-beam  mode, 
the  interferometer  suppresses  radiation  that  is  common  to  both  beams,  regard¬ 
less  of  its  spectrum.  This  includes  diffuse  background  radiation  in  the 
vicinity  of  the  source  of  interest,  absorbing  or  emitting  contaminants  along 
the  line  of  sight  to  the  source,  and  background  (other  than  detector  noise) 
originating  from  within  the  interferometer.  Suppression  ratios  of  100/1  have 
been  achieved  with  our  laboratory  interferometer,  and  there  is  no  obvious 
reason  why  this  cannot  be  readily  improved  to  1000/1  by  future  modifications.  We 
have  also  shown  that  background  structure  with  higher  spatial  frequency  compo¬ 
nents,  that  is  not  eliminated  by  spatial  filtering,  can  be  suppressed  if  its 
spectral  composition  is  distinguishable  from  that  of  the  source  of  interest. 

The  background  suppression  capabilities  of  the  Visidyne  double-beam 
interferometer  with  a  tailored  MTF  have  at  least  two  potentially  important 
applications : 

•  The  remote  sensing  of  effluents,  e.g.,  pollution  monitoring. 

•  Surveillance  and  tracking  systems,  e.g.,  the  detection  of  jet  plumes 
or  missile  trails. 

The  dual  beam  system  is  an  imaging  system  that  can  be  used  with  a  multi¬ 
element  focal  plane  detector  and  it  is  a  diffraction  limited  system  whic^ 


preserves  spatial,  spectral  and  temporal  information.  Two  characteristics  of 
the  double-beam  interferometer  with  tailored  MTF  that  make  it  superior  to 
conventional  detection  systems  include: 

•  Selectable  Spatial-Frequency  Filter 

i)  the  system  spati al -frequency  bandpass  can  be  tailored  to  respond 
to  target  characteristics  while  suppressing  higher  and  lower 
spatial  frequencies. 

•  Real  Time  Optical  Background  Suppression 

i)  reduces  dynamic  range  and  electronic  complexity  required  for 
processing  of  weak  target  signals, 

ii)  reduces  the  effects  of  temporal  variations  in  the  source  or  the 
detectors,  and 

iii)  facilitates  the  extraction  of  a  target  signal  from  background  in 
the  case  of  non-uniform  retroreflecting  mirror  drive 

iv)  suppresses  the  contribution  of  channel  spectra  to  the  measured 
interferogram. 

The  double-beam  interferometer  is  particularly  well  suited  to  the  remote 
sensing  of  stack  effluents.  For  any  given  observation  distance,  the  MTF  of  the 
interferometer  can  be  tailored  to  respond  to  the  volume  of  space  immediately 
surrounding  the  smokestack  while  suppressing  background  radiation  that 
corresponds  to  significantly  larger  or  smaller  angular  sizes.  In  practice, 
regions  of  concentrated  pollutants  with  linear  sizes  from  a  fraction  of  a 
meter  to  tens  of  meters  could  be  observed  with  our  laboratory  interferometer 
at  any  distance  from  ~  100  m  to  tens  of  kilometers  by  varying  the  fore  optics 
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(e.g.,  the  telescope)  and  the  defocusing  lens.  A  mocaic  focal  plane  detector 
is  not  required  for  this  application. 

The  double-beam  interferometer  is  equally  well  suited  for  use  in  surveil¬ 
lance.  The  spatial  filter  can  be  tailored  to  the  size  of  the  plume  if  the 
target  is  measured  in  emission  or  to  the  size  of  the  aircraft  or  missile  if 
the  target  is  measured  in  obscuration. 
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